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Long-term spectral and timing properties of the soft gamma-ray 
repeater SGR 1833-0832 and detection of extended X-ray emission 
around the radio pulsar PSRB1830-08 



P. Esposito/* G. L. Israel,^ R. Turolla,^'^ F. Mattana,^ A. Tiengo,^ A. Possenti,^ 
S. Zane,^ N. ReaJ M. Burgay/ D. Gotz,^ S. Mereghetti,^ L. Stella,^ M. H. Wieringa,^ 
J. M. Sarkissian,^ T. Enoto/'^ R Romano/^ T. Sakamoto/^ Y. E. Nakagawa/^ 
\Q '. K. Makishima/^'^^ K. Nakazawa/^ H. Nishioka/^ C. Frangois-Martin^^ 

^INAF ~ Osservatorio Astronomico di Cagliari, localita Poggio dei Pini, strada 54, 1-09012 Capoterra, Italy 
^INAF - Osserx'atorio Astronomico di Roma, via Frascati 33, 1-00040 Monteporzio Catone, Italy 
^ Universita di Padova, Dipartimento di Fisica, via F. Marzolo 8, 1-35131 Padova, Italy 

'^Mullard Space Science Lahoratoty, University College London, Holmbury St. Maty, Dorking, Surrey RH5 6NT, UK 

^ Frangois Arago Centre, APC (Universite Paris Diderot, CNRS/1N2P3, CEA/DSM. Obserx'atoire de Paris), 13 rue Watt, F -75205 Paris Cedex 13, France 
^INAF — Istituto di Astrofisica Spaziale e Fisica Cosmica - Milano, via E. Bassini 15, 1-20133 Milano, Italy 
^Institut de Ciencies de I'Espai (CSIC-IEEC), Campus UAB, Facultat de Ciencies, Torre C5-parell, E-08193 Barcelona, Spain 
^AIM (UMR 7158 CEA/DSM-CNRS-Universite Paris Diderot} Irfu/Service d'Astrophysique, Saclay; F-91191 Gif-sur-Yvette Cedex, France 
'^Australia Telescope Natioiml Facility, CSIRO Astronomy and Space Science, P.O. Box 76, Epping 1710, Australia 

^° Kavli Institute for Particle Astrophysics & Cosmology (KIPAC), SLAC/Stanford University, PO Box 20450, MS 29, Stanford, CA 94309. USA 
^^INAF — Istituto di Astrofisica Spaziale e Fisica Cosmica - Palermo, via U. La Malfa 153, 1-90146 Palermo, Italy 
^^NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

^^High Energy Astrophysics Laboratory, Institute of Physical and Chemical Research (RIKEN), Wako, Saitama, 351-0198, Japan 
^'^ Department of Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan 
^ ■ Universite Denis Diderot - Paris 7. 4 place Jussieu, F-75252, Paris Cedex 5, France 

m 
(N 
cn 



Accepted 201 1 May 6. Received 201 1 May 4; in original form 201 1 April 12 



o 



X 



ABSTRACT 

SGR 1833-0832 was discovered on 2010 March 19 thanks to the Swift detection of a short 
hard X-ray burst and follow-up X-ray observations. Since then, it was repeatedly observed 
with Swift, Rossi X-ray Timing Explorer, and XMM-Newton. Using these data, which span 
about 225 days, we studied the long-term spectral and timing characteristics of SGR 1833- 
0832. We found evidence for diffuse emission surrounding SGR 1833-0832, which is most 
likely a halo produced by the scattering of the point source X-ray radiation by dust along the 
line of sight, and we show that the source X-ray spectrum is well described by an absorbed 
blackbody, with temperature kT ^ 1.2 keV and absorbing column iVn = (10.4 ± 0.2) x 10^^ 
cm^^, while different or more complex models are disfavoured. The source persistent X-ray 
emission remained fairly constant at ~^3.7 x 10"^^ erg cm~^ s~^ for the first '^20 days after 
the onset of the bursting episode, then it faded by a factor ^40 in the subsequent ^140 days, 
following a power-law trend with index a ~ —0.5. We obtained a phase-coherent timing 
solution with the longest baseline (^225 days) to date for this source which, besides period 
P = 7.565 408 4(4) s and period derivative P = 3.5(3) x 10"^^ s s"\ includes higher 
order period derivatives. We also report on our search of the counterpart to the SGR at radio 
frequencies using the Australia Telescope Compact Array and the Parkes radio telescope. No 
evidence for radio emission was found, down to flux densities of 0.9 mJy (at 1.5 GHz) and 
0.09 mJy (at 1.4 GHz) for the continuum and pulsed emissions, respectively, consistently 
with other observations at different epochs. Finally, the analysis of the field of PSRB1830- 
08 (J 1833-0827), which was serendipitously imaged by the XMM-Newton observations, led 
to the discovery of the X-ray pulsar wind nebula generated by this 85-ms radio pulsar. We 
discuss its possible association with the unidentified TeV source HESS J1834-087. 

Key words: gamma-rays: individual: HESS J 1834-087 - pulsars: general - stars: neutron - 
X-rays: individual: PSRB 1830-08 (J 1833-0827), SGR 1833-0832. 
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1 INTRODUCTION 

Soft gamma-ray repeaters (SGRs; seven confirmed members) and 
anomalous X-ray pulsars (AXPs; twelve confirmed members are 
two classes of X-ray pulsating sources with no evidence for com- 
panion stars which share a number of properties. These include 
rotation periods of several seconds (P ~ 2-12 s), rapid spin 
down (P ~ 10~^^ s s~^), large and variable X-ray luminosi- 
ties (exceeding the rate of rotational energy loss), and the emis- 
sion of flares an d short bursts (see IWoods & Thomps^ l2006t 
iMereghettil boOSi for reviews). SGRs and AXPs (note that the 
distinction is becoming increasingly blurred) are currently inter- 
preted as observational manifestations of magnetars, n amely neu- 
tron st ars powered by their huge magnetic field (e.g. [Paczvnsldl 
19921 : iDuncan & ThompsonI 1 1 992: 'Thomps on & DuncanI 1 1991 



Table 1. Journal of the Swift and XMM-Newton observations. The time of 
the BAT trigger is MJD 55274.774. 



199€ : [Thompson. Lvutikov & Kulkarni 2002). This picture is sup- 



ported by the fact that the dipole magnetic fields inferred for SGRs 
and AXPs from their period and period derivative^ are above, or 
at the high end of, those of the radio pulsars. Surface magnetic 
fields in SGRs/ AXPs in fact often exceed 10^^ G, although an up- 
per limit as l ow as 7.5 x 10^^ G has been recently reported for SGR 
0418+5729 ( lReaetal.ll2010h . 

Between periods of activity, characterised by bursts and signif- 
icant variability in flux, spectrum, pulse shape, and spin-down rate, 
magnetars go through long stretches of quie scence. The discove ry 
of the first 'transient' AXP XTE J1810-197 jlbrahim et al.l2004l) R 
showed that during quiescence magnetars can be faint and not dis- 
similar from hundreds of unidentified sources present in various 
X-ray catalogues (such as the ROSAT, XMM-Newton, and Chan- 
dra ones). This suggested that a potentially large number of mag- 
netars had not been discovered yet and may manifest themselves 
in the future. Thanks also to the effectiveness of the Swift and 
Fermi satellites in catching magnetar outbursts, five new magnetars 
(all of them transients) were discovered in the last few years and 
several major outburs ts were observed from known sources (see 
iRea & Espositd[201 ll and references therein). 

A recent addition to the magnetar family is SGR 1833-0832. It 
was discovered on 2010 March 19 when, at 18:34:50 UT, the Swift 
Burst Alert Telescope (BAT) triggered on a short (<1 s) hard X- 
ray burst and localised it in a region close to the Galactic plan e 
dGelbordetallbOlOl : iBarthelmv et"Zll2010t iGogus et aklboiOah . 
The Swift X-ray Telescope (XRT) started observing the BAT field 
one minute after the trigger and unveiled the existence of a pre- 
viously unknown bright X-ray source. Given the proximity to the 
Galactic plane and the burst properties, the X-ray source was im- 
mediately suggested to be an SGR. The SGR nature of the source, 
now catalogued as SGR 1833-0832, has been confirmed shortly af- 
ter by the discovery with Swif t and Rossi X-ray T i ming Explorer 
(RXTE) of pulsations at 7.57 s l lGogus et al.]|2010bl : lEsposito et al] 
l20ld : IPalmer & GelbordlboiOt) . Subsequent observations also al- 
lowed the determination of the spin-down rate of SGR 1833-0832 
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" This observation was caiTied out in WT mode. 



tesposito et allbOlObl : IGoSus et alj|2010al) . Following the onset of 
the o utburst, the source fl ux remained fairly constant for about 20 
days l lGogtis et alj|201()a ). No radio, optical or infrared counter- 



parts have been det ected Burgav et al.ll201ol : IWieringa et al.ll201(]| : 
iGogus et al.ll2010al) . 

Here we report on the spatial and long-term spectral and tem- 
poral behaviour of the persistent X-ray emission of SGR 1833- 
0832 using new XMM-Newton, RXTE and Swift observations. In 
Sections |2] and [3] we describe the X-ray observations used in 
our study and we present the results of our analysis. In Sec- 
ti on [4] we give more detai ls on the radio observations presented 
in lBtjrgav et al] ilOld ) and lWieringa et all ( l2010h . In Section[5]we 
report on the analysis of the field of PSRB 1830-08 (J 1833-0827), 
which was serendipitously imaged by the XMM-Newton observa- 
tions, and discuss its possible association with the unidentified TeV 
source HESS J1834-087. Discussion follows in Section|6] 



See the McGill Pulsar Group catalogue at the webpage 
http://www.physics.mcgill.ca/~pulsar/magnetar/main.html. 
^ According to the usual magnetic braking formula, for a neutron star of 
10 km radius and 1.4 sol ar masses, Bjip ~ 3.2 X 10l^(PP)^/2 G (e.g. 
iLorimer & Krameil2004l) . 

Actually, at that time another magnetar with a transie nt behaviour, 
SGR 1627^1 teouveliotou et alJl2003l : lEsposito et alj|2008l) . was already 
known, but its chai'acteristics were not well established yet. Signiticant 
flux variability was reported also in the AXP candida te AX J1845. 0-0300 
iXorii et aljigillVasisht et alj2000tlTam et al.ll2006l) . 



2 X-RAY OBSERVATIONS 
2.1 XMM-Newton 

The three foc al plane CCD came ras of the XMM-New ton EPIC in- 
strument, pn jStriider et allbOOll) . MOS 1 and M0S2 jTumer et al] 
l200lh . cover the 0.1-12 keV energy range with an effective area 
of roughly 1400 cm^ for the pn and 600 cm^ for each MOS. After 
the discovery of SGR 1833-0832, XMM-Newton pointed its mir- 
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rors towards the new SGR three times (see Table [T). All observa- 
tions were performed with the thick optical filter and in full frame 
modeQ except for the first one, carried out with the MOS cameras 
in large window mode. Moreover, XMM-Newton serendipitously 
imaged the field of SGR 1833-0832 on 2006 September 16 during 
an observation targeting the nearby supernova remnant G23. 5-0.0 
(obs. ID: 0400910101, exposure: 12.4 ks; all the detectors were in 
full frame mode with the medium filter). 

The data were processed using version 10.0 of the XMM- 
Newton Science Analysis Software (SAS) and standard screening 
criteria were applied. Source events were accumulated for each 
camera from circular regions with a 36 arcsec radius. We selected 
this aperture, corresponding to ~85% of the encircled energy frac- 
tion at 5 keV for a point source, in order to minimise the contami- 
nation from the diffuse emission surrounding SGR 1833-0832 (see 
Section l3Tt . The background counts were extracted from source- 
free regions far from the position of the SGR. The ancillary re- 
sponse files and the spectral redistribution matrices for the spectral 
analysis were generated with the SAS tasks ARFGEN and RMFGEN, 
respectively. 

2.2 Swift 

The X-Ray Telescope (XRT: [Burrows et^ IHool) on-board Swift 
uses a CCD detector sensitive to photons between 0.2 and 10 keV 
with an effective area of about 1 10 cm'^. Twenty-seven observations 
of SGR 1833-0832 were performed (see Table [T), starting right af- 
ter its discovery, in b oth photon count ing (PC) and windowed tim- 
ing (WT) modes (see lHill et alll 20041 for more details on the XRT 
readout modes). The res ults of the firs t 18 Sw ift observations have 
already been published in lG6gusetal.l ( l2010al) . while the others are 
reported here for the first time. 

The data were processed and filtered with standard criteria 
using the FTOOLS software package. We extracted the PC source 
events from a circle with a radius of 15 pixels (one pixel corre- 
sponds to about 2'.' 36) and the WT data from a 25-pixel- wide strip. 
To estimate the background, we extracted PC and WT events from 
source-free regions distant from the position of SGR 1833-0832. 
For the spectral fitting we used the latest available spectral redistri- 
bution matrix in CALDB (vOll), while the ancillary response files 
were generated with XRTMKARF, and they account for different ex- 
traction regions, vignetting and point-spread function corrections. 

2.3 RXTE 

The Proportional Counter Amy (PCA; Ijahoda et all 1 19961) on- 
board RXTE consists of an array of five coUimated xenon/methane 
multi-anode Proportional Counter Units (PCUs) operating in the 2- 
60 keV energy range, with a total effective area of approximately 
6500 cm'^ and a full width at half-maximum field of view (FOV) of 
about r. The 80 RXTE/PCA pointed observations of SGR 1833- 
0832 (obs. ID: 95048) reported here span from 2010 March 19 to 
2010 December 04. The exposures range from 1 ks to 20 ks, for 
a total of about 466.4 ks. The results of the first 3 1 RXTE obser- 
vations have already been published in iGogus et al.l ( l2010ah . while 
the analysis of the others is reported here for the first time. 

The raw data were reduced using the FTOOLS package. Given 
the non-imaging nature and wide FOV of the PCA instrument and 

See the XMM-Newton User Handbook at 
|http://xmm.esac.esa.int/extemal/xmm_user_support/documentation/index.shtml| 
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Figure 1. Surface brightness radial profiles of the background-subtracted 
X-ray emission of SGR 1833-0832 from the XMM-Newton MOSl and 
MOS2 cameras (3-6 keV data from the three observations combined). 
The red solid line sh ows for each camera the best-fit point-spread func- 
tion jGhizzardil200lh . The radial profiles show clear evidence for extended 
emission starling from ^35 ai'csec. 



the relatively low flux of SGR 1833-0832, the RXTE data were used 
only to study the timing properties of the source. We thus restricted 
our analysis to the data in Good Xenon mode, with a time resolution 
of 1 fis and 256 energy bins. They were extracted in the 2-10 keV 
energy range from all active PCUs (in a given observation) and all 
layers, and binned into light curves of 10-ms resolution. 



3 X-RAY DATA ANALYSIS AND RESULTS 

We inspected all observations for the presence of bursts by a care- 
ful examination of the light curves binned with differ ent time res- 
olution s. A few were found (besides those reported in lGogiis et al.l 
l2010al we found another short and weak burst in the RXTE data of 
2010 May 18). In the analyses that follow, we removed the bursts 
from the event lists by applying intensity filters. Swift and RXTE 
data from a few contiguous observations carried out with the same 
instrumental setup were combined in order to achieve better statis- 
tics and higher signal to noise ratio. For the timing analysis, the 
data were corrected to the barycentre of the Solar system u sing the 
ChandmfUKlKT position reported in iGogus etalj|2010al) . 



3.1 Spatial analysis 

As can be seen in Fig.[T] some diffuse emission around SGR 1833- 
0832 was detected during the three XMM-Newton observations, 
mainly in the 3-6 keV energy band. Considering the large ab- 
sorption derived from the X-ray spectrum, this is likely a halo 



produced by the scattering o f the point source X-ray radiation by 
dust along the line of sight llOverbecklll96^ : iTiengo et al.l 120101 : 
iRivera-Ingraham & van Kerkwiikll201(t) . 

By analysing the spectrum from different annular regions, we 
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have found that, as expected from a dust halo, the diffuse X-ray 
emission is significantly softer than the point source. Moreover, the 
XMM-Newton observation of the SGR 1833-0832 field performed 
in 2006, where no point-like nor diffuse emission was detected at 
the SGR position (see Section [J!4t . favours the dust halo hypothesis 
with respect to steady diffuse emission, as expected, for example, 
in a supernova remnant. However, the data quality of the outburst 
observations is not good enough to detect significant variability in 
the extended emission at the level expected for a dust scattering 
halo. 



3.2 Spectral analysis 

For the spectrosco py (performed with the XSPEC 12.6 fitting pack- 
age; |Am^2[l99^ we concentrate first on the XMM-Newton spec- 
tra which, owing to the EPIC instrument high throughput and long 
observing time, are those with the best statistical quality. We fit 
the spectra from the three observations simultaneously with the hy- 
drogen column density tied between all data sets. Photons having 
energies below 2 keV and above 10 keV were ignored, owing to 
the very few counts from SGR 1833-0832. This resulted in about 
6100 ± 80, 5800 ± 90, and 4000 ± 70 net EPIC-pn counts in the 
three obs ervations in chronological o rder. The abundances used are 
those of I Anders & G revesse' ("1989^ and photoelectric absorptio n 
cross-sections are from Balucinska- Church & McCammonI l ll992h . 

Using the first twelve 5w(/f/XRT observations in PC mode, 
iGogiis et al.l ( 1201 Oal) could not discriminate between a power-law 
and a blackbody as the model that better describes the spectrum 
of the SGR persistent emission. A fit of the XMM-Newton pn and 
MOS data with an absorbed power law yields a relatively high 
value [xl = 1.21 for 673 de grees of freedom (dof)] and structured 
residuals, while a better fit (xS = 1-05 for 673 dof) is obtained 
using a blackbody model (corrected for the absorption). No addi- 
tional spectral components are statistically required. The results of 
this simultaneous modelling are presented in Table|2](the values of 
the spectral parameters were not significantly different when each 
observation was fit separately). The best-fitting hydrogen column 
density is (10.4 ± 0.2) x 10^^ cm~^ (here and in the following 
uncertainties are at Icr confidence level, unless otherwise noted) 
and the blackbody temperature is consistent within the errors with 
being constant at kT ~ 1.2 keV, while the observed flux slightly 
(but significantly, see Table© decreased from F ~ 3.9 x 10"^^ 
erg cm"^ s"^ to 3.1 x 10"^^ erg cm"^ s"^ (in the 2-10 keV band) 
during the ^-^20 days that separate the first and last XMM-Newton 
pointings. Assuming an arbitrary distance of 10 kpc, the corre- 
sponding XMM-Newton maximum and minimum luminosities are 
~9x 10^* erg s'^ and ~7x 10^* erg s'^ (Tabled. For each XMM- 
Newton observation we performed phase-resolved spectroscopy by 
extracting the spectra for different selections of phase intervals. No 
significant variations with phase were detected, all the spectra be- 
ing consistent with the model and parameters of the phase-averaged 
spectrum, simply re-scaled in normalisation. 

Most magnetars, especially when in outburst, exhibit more 
complex spectra, usually fit by the superposition of two/three 
blackbodies or a blackbody and a high- energy p ower law (e.g . 
Mereghetti et al.l2005l;[Esposito et al.l2()09 ; Bemar dini et al.l2009l 



0.01 




20111 ; llsraeletal.1 l2007al ; lEnoto et a'nT2009t iRea et al.l I2OO9I) . 
The spect rum of SGR 1 833-08 32 is reminiscent of that re- 
ported by lEsposito et al.l ( l2010ah for SGR 04 18-1-5729 during its 
2009 outburst . Howev er, at the time of the analysis reported in 
lEsposito et alj ( 1201 Oal) . only low-counts-statistics Swift and RXTE 
spectra were available for SGR0418-F5729. Higher quality XMM- 



2 5 10 

Energy (keV) 

Figure 2. Fit of the spectrum obtained from the first two XMM-Newton 
obsei'vations (see Section 13.21 with the single-blackbody model. Black 
squares, red circles, and blue stars represent the pn, MOS 1 , and MOS2 data, 
respectively. Bottom panel: the residuals of the fit (in units of standard de- 
viations). 



Newton data, now public, clearly show that a more complex model 
is required for the emission of SGR0418-I-5729I!] 

In order to achieve better statistics and higher signal to noise, 
we merged the data from the first two XMM-Newton observations 
(since there is no evidence for variations in the flux and spectrum of 
the source between them) and accumulated a combined spectrum 
which is presented in Fig. |2] Again, the single-blackbody model 
provides an excellent fit (xt = 1-00 for 352 dof), and no additional 
components are required. By including in the m odel a power-law 
component with photon index fixed to 3 (see e.g. lRea et al.ll200^ . 
we can set 3a upper limits of ~25% and «30% on the contribu- 
tion of this component to the total observed and unabsorbed fluxes, 
respectively. This shows that the relatively low flux and high ab- 
sorption of SGR 1833-0832 make even our deep XMM-Newton ob- 
servations not very sensitive to the presence of a second spectral 
component in this SGR. We repeated the phase-resolved spectral 
analysis on the combined dataset. Once more, we see no evidence 
of spectral shape evolution with the rotational phase; in particular, 
no variation larger than the la error (kST ~ 0.02 keV) was found 
in the blackbody temperatures at different phases. 

To obtain flux measurements over the outburst, the 5wi/f/XRT 
data were fit in the same way as the XMM-Newton ones, simultane- 
ously and with all parameters left free to vary, except for the absorp- 
tion column density, that this time was fixed at the value measured 
with XMM-Newton. This resulted in an acceptable fit (xl = 1-01 
for 602 dof) with spectral parameters similar to those reported in 
Table|2](see Fig.[3]l. We plot the resulting long-term light curve in 
the top panel of Fig. [3] 

Although the available data do not allow us to perform an ac- 
curate modelling of the decay shape, because of the relatively high 



We have analysed an XMM-Newton observation of SGR 041 8+5729 per- 
formed on 2009 August 12-13, about 68 days after the outburst onset (obs. 
ID 0610000601, exposure 67.2 ks). The data are well described by a power 
law plus blackbody model (xl = 1-1 for 231 dof) with the following 
spectral parameters: A'^jj = (6.4 ± 0.3) X 10^^ cm~^, photon index 
r = 2.6lo^, blackbody temperature kT = 0.93 ± 0.01 keV, and 0.5-10 
keV absorbed flux (7.0 ± 0.2) X IQ-^^ gj-g ^.^-2 ^-1 ^ double black- 
body gave a poorer fit (x^ = 1.2 for 231 dof), whereas both blackbody and 
power-law models yield statistically unacceptable fits (x^ ^ 2). 
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Table 2. Spectral fit of the XMM-Newton data with the blackbody model (xj = 1.05 for 673 
dof). Errors are at a Ic confidence level for a single parameter of interest. The measured Afjj 
value is (10.4 ± 0.2) X 10^2 cm^^ 



Obs.ID 


kT 
(keV) 


Radius" 
(km) 


Absorbed flux'' 
(10^12 ergcm-2 s"!) 


Luminosity'^ 
(1034 ergs-i) 


0605851901 
0605852001 
0605852101 


1.20 ±0.01 


0.62 ±0.02 
0.62 ±0.02 
0.55 ±0.02 


3.9 ±0.1 
3.6 ±0.1 
3.1 ±0.1 


8.91°;? 

8.4 ±0.1 
7.2 ± 0.1 



" The blackbody radius is calculated at infinity and for an arbitrary distance of 10 kpc. 
^ In the 2-10 keV energy range. 

In the 2-10 keV energy range and for an arbitrary distance of 1 kpc. 



uncertainties on the fluxes and moderate time-span, we observe that 
similarly good fits can be obtained with either an exponential func- 
tion or a broken power-law model. For an exponential function of 
the form F{t) = Aexp{-t/r), the best-fitting values {xl = 0.84 
for 23 dof) are = (3.92 ± 0.05) x 10"^^ ergcm"^ s"^ and 
T = (108 ± 9) day0 Adopting a broken power law (xi ~ 0.73 for 
22 dof;), the break occurs at (20 ± 3) d, when the index changes 
from Qfi = -0.01 ± 0.02 to 02 = -0.54 ± 0.09; the flux at the 
break time is (3.7 ± 0.2) x 10~^^ erg cm^^ s~^. In both cases we 
assumed as t = the time of the 5w//f/BAT trigger. From the bro- 
ken power-law fit in particular, it is apparent that the source flux 
is consistent with a constant value for the first ~20 days. The two 
models considered are plotted in Fig. |3] 



3.3 Timing analysis 

In lEsDosito et al.U2010l) we reported on our search for periodicities 
made on the first SwiftlXKY observation (00416485000) by calcu- 
lating a fast-Fourier-transform power spectrum. A very prominent 
peak occurs in the spectrum of that observation at 7.5653(4) s (the 
quoted uncertainty indicates the Fourier period resolution). Pulsa- 
tions were clearly detected also in all the others Swift and XMM- 
Newton datasets, and in the RXTE ones up to 2010 October 30 
(MJD 55499), when presumably the flux became too low for the 
PCA sensitivity. 

In order to obtain a refined ephemeris for the longest pos- 
sible baseline, we studied the pulse phase evolution in these 
observati ons by means of an iterative phase-fitting technique 
(see e.g. iDall'Osso et al. II200I). The fits were carried out in the 
range 2-10 keV with a minimisation approach using MINUIT 
jjames &Roo3ll975h . 

Throughout the period covered by useful observations (~225 
days) the relative phases and amplitudes were such that the phase 
evolution of the signal could be followed unambiguously. A 
second-orde r polyno mial, as employed in the recent analysis by 
iGogiis et al.l ( 1201 Oal) over the first ~47 days, provides an unac- 
ceptable fit to the data, with xt = 9-06 for 86 dof. We tried 
higher-order polynomials until the addition of a further (higher- 
order) term was not statistically significant at more than 3 a with 
respect to the null hypothesis (as evaluated by the Fisher test). 
The outcome of this process was a fourth-order polynomial (the 
improvement obtained in the fit with a fifth-order polynomial has 



^ An additional component Fq representing the base flux level [F{t) = 
Fq + Aexp(— f/r)] is not required. Including it anyway in the fit (xj — 
0.83 for 22 dof), we find Fq = (5 ± 4) X 10~" erg cm-^ j-i. 




Time (days since MJD 55274.77418981) 



Figure 3. Time-evolution of the characteristic parameters of SGR 1833- 
0832 inferred from the spectral analysis (single blackbody model) of the 
Swift (black triangles) and XMM-Newton (red circles) data. The different 
panels show (from top to bottom): the (absorbed) flux in the 2-10 keV en- 
ergy range (two possible models describing the decay are also plotted: an 
exponential function with the solid magenta line and a broken power-law 
with the dot-dashed blue line; see Section 13. 2t . the blackbody tempera- 
ture and radius (evaluated at infinity and assuming an arbitrary distance 
of 10 kpc), and the pulsed fraction. We assumed as i = the time of the 
Swift/BAT trigger. 



a statistical significance of only 2.6o-), which we used to fit the 
phase shifts. The resulting phase-coherent solution is given in Ta- 
ble|3]and plotted in Fig.lH the best fit (xl = 1-07 for 84 dof) gives 
u = 0.132180 571(7) Hz and i> = -6.0(5) x 10"" Hz s"\ 
assuming MJD 55274.0 as ref erence epoch. We ha ve checked that 
the positional uncertainty (0V3; lGogtis et al.l2010ah does not signif- 
icantly affect the rotational parameters resulting from our analysis. 
In Fig. |5] we show the three XMM-Newton light curves ob- 
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Table 3. Spin ephemeris of SGR 1833-0832 obtained using the combined 
RXTE, Swift, and XMM-Newton observations. We also give for conve- 
nience the coiTesponding peiiod P and period derivative P, as well as 
the derived characteristic age Tc = P/{2P), dipolar magnetic field B 
(Sc^/PP/Svr^Rf^) 1/2, and rotational energy loss E = iw^IPP-'-' (here 



we took = 10 km and / : 
of inertia, respectively). 



10'*^ g cm'^ for the star radius and moment 



Parameter 


Value 


Range (MJD) 


55274.775-55499.170 


Epoch (MJD) 


55274.0 




0.132180571(7) 


V (Hz s-i) 


-6.0(5) X lO-i* 


i>(Hzs-2) 


-1.3(2) X 10^2" 


V (Hz s-3) 


9(2) X 10-28 


X^/dof 


89.47/84 


P(s) 


7.565 408 4(4) 


P(s s-i) 


3.5(3) X 10-12 


Tc (kyr) 


35 


B(G) 


1.6 X 10" 


£:(ergs-i) 


3.2 X 10^2 




50 100 150 

Time (days since MJD 55274.0) 



200 



Figure 4. Top panel: Swift, XMM-Newton, and RXTE pulse phase evolution 
with time with respect to the period measured during the first Swift observa- 
tion. The solid red line represents the best-fitting fourth-order polynomial 
solution. Bottom panel: Time residuals with respect to timing solution. 



tained folding the high time-resolution (nominal frame time of 
73.4 ms) EPIC-pn data at our phase-coherent ephemeris. The pulse 
profile is sinusoidal and the pulsed fraction, that we define as 
the semi-amplitude of sinusoidal modulation divided by the mean 
source count rate, was consistent through the first two observa- 
tions [(69.0 ± 1.4)% and (67.4 ± 1.5)%, respectively], while it 
decreased to (57.0 ± 1.7)% in the third one. We have investigated 
the morphology of the pulse phase distribution as a function of en- 
ergy by comparing the XMM-Newton pulse profiles in different en- 
ergy bands. The measured pulsed fractions in the three observations 
are (66 ± 2)%, (64 ± 2)%, and (56 ± 2)% in the soft (2-5 keV) 
energy band, and (75 ± 2)%, (71 ± 2)%, and (59 ± 3)% in the 
hard (5-10 keV) band. Apart from this marginal indication for an 
increasing trend of the pulsed fraction with energy, no significant 
pulse shape variations (such as phase shifts of the maxima) were 
found as a function of energy by cross-coiTelating or comparing 
through a two-sided Kolmogorov-Smimov test the soft and hard 
folded profiles. 

The 2-10 keV pulsed fractions measured in the individual ob- 
servations obtained with XMM-Newton and Swift are shown in the 
bottom panel of Fig. [3] (we did not consider the RXTE data as the 
non-imaging PCA instrument does not ensure reliable background 
subtraction). A constant fit of the pulsed fraction values derived 
with XRlQ (and shown in Fig. [3} does not adequately describe the 
data (xt = 4.52 for 24 dof); however no particular trend is ap- 
parent in the pulsed fraction evolution and the simple functions we 
tried do not yield significantly better fits. 



counts (we did not use the MOS data since the source position 
fell near CCD gaps in both cameras). The net exposure time was 
9.2 ks. Considering the expected background counts (estimated 
from large suiTounding source-free regions) and assuming Poisso- 
nian fluctuations, we set a 3cr upper limit on the 2-10 keV count- 
rate of SGR 1833-0832 of 2.3 x 10"^ counts s"\ Assuming a 



blackbody spectrum with kT 

-,22 „„-2 



1.2 keV and an absorption of 



the 2-10 keV observed flux of ~4 x 10" 



-14 —2 - 

erg cm s 



and of 



~8 X 10~ erg cm- s^ on the 2-10 keV unabsorbed flux. If in- 
stead we assume a much softer blackbody spectrum with kT = 0.3 
keV, as it might be expected from a magnetar in quiescence (e.g. 
iLamb et aljr2002l : lOotthelf et al.ll2004 ISemardini et al.ll2009l) , the 
following upper limits are derived (2-10 keV band): ~2 x 10~^^ 
erg cm-2 s^^ on the observed and ~2 x lO"'^^ erg cm-^ on 
the unabsorbed flux. 

More recently, on 2009 February 13, the field of SGR 1833- 
0832 wa s serendipitously imaged for 8 ks by th e Chandra/ ACIS-S 
(see alsolMisanovic. Kargaltsev & Pavlovl201 ih . The SGR was not 
detected and Gog u fj, et al. (2010a|) ' derived a 2a upper limit on the 
absorbed flux of 3.4 x 10"^'^ erg cm-^ s'^, ass uming a power-law 
mode l with F - 3 and iVn 14 x 10^^ cm-^ jGelbord & Veterd 
I2OIOI) . 



3.4 Archival searches 

No X-ray sourc e was detected at the Chandra position of 
SGR 1833-0832 jCogus et al.| |2010a| ) in the XMM-Newton obser- 
vation performed in 2006. In a 20-arcsec circle centred at the SGR 
coordinates we detected in the 2-10 keV energy band 24 EPIC/pn 



^ We did not attempt to fit simultaneously the XMM-Newton and Swift data 
in order to avoid possible effects due to the different responses and energy 
dependence of the effective areas of the of the EPIC and XRT detectors. 



4 RADIO OBSERVATIONS 
4.1 ATCA 

We observed the Swift/XKT position of SGR 1833-0832 at 20 cm 
twice, on 2010 March 22 (UT time 21:40-22:20) and March 23 
(UT time 00:17-01:32), with a total on-source integration time of 
105 minutes. The observing band was 1.1-1.9 GHz, split in 1 MHz 
channels and with a central frequency 1.52 GHz; however, several 
channels, affected by radio frequency interferences, were excised 
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Figure S. Background-subtracted XMM-Newton epoch-folded pulse pro- 
files (32-bin, EPIC-pn). The sinusoidal fit to the data is superimposed in 
red (the values are 0.9, 1 .0, and 1.2 for 29 dof, for the three observations 
in order). 



from the spectrum. The array was in the hybrid configuration H 168 
(maximum baseline ~4. 5km)|f| 

Given the location of the source in the Galactic plane, images 
including the short baselines were too confused by the emission of 
extended sources to yield useful limits. Continuum images combin- 
ing all usable frequency channels were then made considering only 
the long baselines, i.e. those including the antenna CA06. These 
images were still somewhat confused by other sources in the field 
and the rms noise near the field centre was twice that outside the 
X 80" wide beam. In particular, the local rms at the SGR 
position turned out to be 0.3 mjy, whence the 3ct flux density limit 
for an unresolved continuum radio emission from SGR 1833-0832 
is 0.9 mJy. 



4.2 Parkes 

The source was observed on 2010 March 25 (MJD 55279.954) 
for 90 minutes at a centre frequency of 1374 MHz with the 
central beam of the Parke s 20 cm multibeam (MB) receiver 
dStavelev-Smith et all [199^ . The total 288 MHz bandwidth was 
split into 96 channels, each 3 MHz wide, and the time series was 1 
bit sampled every 1 ms. Data have been searched for periodicities 
around the value obtained from the X-ray observations over a wide 
range of dispersion measure values (DM < 6000 pc cm^'^, given 



the very high A^'h value derived from X-ray observations). A search 
for single de-dispersed pulses was also carried out. 

No radio pulsation with a period matching (within ±2.5 ms) 
the X-ray period, nor half of it, has been found. The single pulse 
search only revealed the presence of a bright known pulsar in the 
beam (the 85-ms J1833-0827/B1830-08, at DM 411 pc cm"^'; see 
Section Is} ■ The upper limit on the pulsed flux density at 20 cm 
wavelength is of 0.09 mJy for a duty cycle of 10 per cent, calcu- 
lated for the maximum DM investigated and assuming a minimum 
signal-to-noise ratio of 9. The latter relatively high value for the 
adopted minimum signal-to-noise ratio is due to the presence of 
strong radio interferences at ~7.5 s, that affect the signal folded at 
the SGR period. 



5 X-RAY EMISSION AROUND PSR B1830-08 

Inspection of the XMM-Newton EPIC images reveals an excess of 
X-ray emission north of SGR 1833-0832, encompassing the po- 
sition of the midd le-aged radio pulsar PSR B 1830-0 8 (spin-down 
age Tc ~ 150 kvr: lTavlor. Manchester & Lvnelll993l) . 

PSR B 1 830-08 is located outside, and possibly associated to, 
the asymmetric shell-type supernova remnant W41 (G23. 3-0.3) 
despi te an angular separation of 24' to the ce ntre of the rem- 
nant ( lKassinjfl993 : iGaensler & JohnstorJll995l) . Its distance es- 
timated through the dispersion measure (DM — 411 pc cm"^) 
and the NE2001 model of the Galactic electron density distribu- 
tion bv ICordes & Laziol ( l2002h is 4.6 kpc. PSR B 1830-08 might 
be also linked to the unidentified TeV sour ce HESS J 1834 -087 ly- 
ing inside W41 ( Aharonian et al. 2006; Al bert et al.ll2006l) . In this 
case, the extended TeV source would be a relic pulsar wind nebula 
(PWN) generated by electrons ejected by PSR B 1 830-08 at its birth 
place and up-scattering th e surrounding low-energy photons (e.g., 
iBartko & Bednarekl l2008h . However, different explanations have 
been suggested: the TeV emission of HESS J 1834-087 could result 
from the tt" decay of shock-accelerated hadrons interacting with 
a giant molecular cloud along the line of sight jTian et ajll2007h 
or from the inverse Compton radiation of electr ons in a different 
PWN rece ntly discovered near the centre of W 41 dMukheriee et al.l 
12009; M isanovic, Kargaltsev & Pavlo^^l201 l l ) . The detection of a 
large-scale X-ray nebula trailing the pulsar in the direction of the 
supernova remnant and of the TeV source would support the asso- 
ciation between PSRB1830-08 and HESS JI834-087. 

We refined the imaging analysis by considering only EPIC- 
MOS data from the observations 0605852001 and 0605852101, 
whereas PSR B 1 830-08 was located in the unread regions of the 
MOS detectors operating in Large Window mode during obser- 
vation 0605851901. The EPIC-pn data of the three observations 
were not used either, because of a higher level of residual soft 
proton contamination indicated by the ratio of the surface bright- 
ness inside the field of vie w to the one outside the field of view 
toe Luca & Molendill2004l) . We selected only 1- and 2-pixel events 
with the default flag mask. 

We built background maps using the Filter Wheel Closed 
(FWC) event lists provided by the XMM-Newton background work- 
ing grouplfl As the profile of the FWC maps changed since the 
launch of the mission, we filtered only the events collected after 
Januaiy 2007 to produce a more accurate template of the particle- 
induced background. The background maps were re-scaled in nar- 



* See |http://vi'Vi'Vi'.narrabri.atnf.csiro.aii/observing/configs.html| for more 
details. 



See |http://xmm.vilspa.esa.es/extemal/xmm_svi'_cal/background/index.shtml| 
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Figure 6. Left: XMM-Newton EPIC-MOS mosaiced image. The image is background-subtracted, exposure-corrected and adaptively smoothed; the intensity 
scale has been chosen to highlight the faint diffuse emission around PSRB1830-08. The brightest source in the centre is SGR 1833-0832. The pulsar position 
is marked with a cross. The arrow indicates the pulsar's measured direction of motion, the length of the arrow corresponding to the distance travelled by the 
pulsar over 10 kyr. The two solid lines from the pulsar position to east indicate the extrapolation of the pulsai' proper motion backwards over the last 50 kyr 
accounting fo r the measure uncertainties. The large solid circle e ncloses the EPIC field of view. The extensions of the supernova remnant W41 (dashed circle; 
lKassim|[l992h and of HESS J 1834-087 (Icr, small solid circle; lAharonian et al.ll200^ are also shown. Right: zoom-in of the EPIC-MOS mosaiced image 
shown in the left panel. The contours are uniformly spaced. The inner and outer dashed ellipses outline the 'compact' and 'halo' emissions, respectively. The 
length of the compass' arms is 2'. 



row energy bands according to the counts in the non-vignetted 
MOS comers. Finally, the images were mosaiced, background- 
subtracted, exposure-corrected, and adaptively smoothed to achieve 
a signal-to-noise ratio of 6. 

Diffuse emission around PSRB 1830-08 is apparent in the 2- 
10 keV band mosaic, extending ~2' north-east and south of the 
pulsar (Fig. |6] left panel). A closer inspection shows a compact 
structure at the pulsar position embedded in a halo of lower sur- 
face brightness (Fig. |6] light panel). The compact structure is ~l' 
elongated towards south-east, whereas the halo is mostly extending 
in the perpendicular direction. The bulk of the halo emission lies 
closer to PSRB 1830-08, while the the south-western tip is appar- 
ently connected with SGR 1833-0832; this is most likely due to the 
smoothing procedure. Both the halo and the compact structures are 
hardly visible in the soft band (0.5-2 keV), suggesting an absorbed 
non- thermal emission. This is confirmed by the spectral analysis. 

We produced the spectra also including the three EPIC-pn data 
sets. The source counts were extracted from the mutually exclusive 
elliptical regions shown in Fig. |6] (right panel). The background 
counts were estimated from 60" apertures placed on the same CCD 
of each camera, and at the same distance to the readout node as the 
source for the EPIC-pn. For each EPIC camera, we merged the in- 
dividual spectra from the different observations into a single spec- 
trum and combined the response matrices. The compact and halo 
regions contained 560 ±30 and 2800 ± 160 background-subtracted 
counts, respectively, corresponding to an average surface bright- 
ness of 0.2 and 0.085 counts arcsec"'^. The spectra were rebinned 
with a minimum of 20 counts per spectral channel for the compact 
emission and 50 for the halo. 

The EPIC pn and MOS spectra of the compact structure (Fig. 
|7J are well fit by an absorbed power-law model, yielding photon 
index F = l.Olg ^, unabsorbed 2-10 keV flux f(2-iokcV) ~ 



i^-^-o t) ^ 10 erg cm ^ s ^, and hydrogen-equivalent column 
density 7Vh= (4 ± 2) x 10^^ cm"^ (reduced = 1-0 for 46 
dof). The halo spectra are also well fit by an absorbed power-law 
model withr — 1.7lo4, unabsorbed 2-10 keV flux F(2_iokcV) = 
(S.glo.g) X 10"^=* ergcm'2s~\ and A^h^ (2.9lJ;g) x 10^^ 
cm~^ (reduced ~ 1-1 for 153 dof). By fitting simultaneously 
the compact and halo spectra with the same column density, we 
found iVH= (3.4l|!,;g) x 10^^ cm~^, slightly higher than the neu- 
tral hy drogen column density fro m radio observations ~2 x 10^^ 
cm l lDickev & LockmarJll99(3l) , which might indicate additional 
absorbing material along the line of sight. 

The hard diffuse emission around PSRB 1830-08 is sugges- 
tive of a PWN. With a spin-down power of ij = 5.8 x 10^^ erg s~ \ 
PSRB 1830-08 is energetic enough to produce a synchrotron neb- 
ula detectable in X-rays for a broad range of distances. For an as- 
sumed distance of 4.6 kpc, the measured fluxes imply a luminosity 
in the 2-10 keV energy range of ~4 x 10''^ erg s~^ for the compact 
structure and of ~1.5 x 10'^'' erg s~^ for the halo. This corresponds 
to a conversion efficiency Lx/E of 7 x 10~^ and 3 x 10~'^, within 
the observed range of pul sar/PWN systems dPossenti et al.ll20o3 ; 
iKargaltsev & Pavlovll2008l) . 

The inner compact stmcture looks one-sided, and its bright- 
ness peak is coincident with the pulsar position. Such a morphol- 
ogy is reminiscent of a bow-shock PWN, in which the wind of a su- 
personically moving pulsa r is confined by the ram- pressure of the 
surrounding medium (see lOaensler & Slanell2006l for a review). 
This phase occurs once the pulsar escapes its remnant, several 
ten thousand years after the supernova event, which is compatible 
with the age of PSR B 1830-08. Furthermore, the proper motion of 
PSRB 1830-08 (33 ± 5 mas yr~ ^ , corresp onding to a linear veloc- 
ity of 730ci4.6 km s' MHobbs et alj|2004) lies in the high-velocity 
end of the pulsar population and it is most likely supersonic. The 
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Figure 7. EPIC MOSl (red circles), MOS2 (blue stars), and pn (black 
squares) spectra of the compact structure around PSR B 1 830-08 (inner el- 
lipse in Fig.|6] right) fit to the model described in the text. The lower panel 
shows residuals from the best fit in units of Icr. 



compact emission seems to be aligned with the pulsar proper mo- 
tion (Fig. [6] right), the brightness profile declining gradually from 
the pulsar position to south-east and sharply to north-west. 

The location opposite to the direction of the pulsar motion, 
the hard extended emission, and the column density suggest that 
the halo is also linked to PSR B 1830-08. However, its morphol- 
ogy is dif ferent from observed (e.g., the 'Mouse' G359.23-0.82 , 
lGaensler"et al. 2004) and simulated (e.g.. iBucciantini et al.ll200^ 
bow-shock PWNe, in which the emission is tightly confined along 
the proper motion direction. The halo forms instead a broad diffuse 
region behind the pulsar, almost perpendicular to the proper motion 
direction. 



6 DISCUSSION 
6.1 SGR 1833-0832 

In this work, we presented a new timing solution for SGR 1833- 
0832 based on a phase-coherent analysis of XMM-Newton, Swift, 
and RXTE data and valid in the range from MJD 55274.775 to 
55499.170. Under the standard assumption that the neutron star 
slows down because of magnetic braking, it implies a character- 
istic age Tc — 35 kyr, a dipolar magnetic field B ~ 1.6 x 10^'* 
G and a rotational energy loss rate E ~ 3.2 x 10^'^ erg s^'^ (see 
Tabled. 

The modelling of the phase shifts required a fourth-order poly- 
nomial which in general is an indication for timing noise ('poly- 
nomial whitening'). However, it is worth noting that the second 
derivative we measured is unlikely to be related to a random change 
of the pulse profiles, which is expected to introduce only a ran- 
dom distribution of the phase residuals, rather then a cubic term. 
Recent studies on a sample of 366 radio pulsars showed that cu- 
bic terms in the phase residuals are actually possible but smaller 
by several orders of magnitudes than those detected in SGR 1833- 
083f°l and recorded on time-scales longer (years) than those we 



^" The indicator commonly used to quantify th e amount of timing noise i n 
radio pulsars is the Ag parameter introduced bv lArzoumanian et al.l ( ll994l) . 
The value computed with the ephemeris of SGR 1 833-0832 given in Ta- 
ble[3]is Ag ~ 4.2, which is much higher than t he typical values measured 
for radio pulsars (roughly between and —5: lArzoumanian et al.|[l994l ; 



are sampling in our dataset jHobbs. Lvne & Kramej2oTol) . We also 
note that the (long-term) timing noise of young radio pulsars (age 
< 10' yr) can be better unde rstood as resul ting from the recovery 
from previous glitch events jHobbs. Lvne & Kramer 2010), mak- 
ing it unlikely that the frequency second derivative is actually due 
to random noise. 

A possibility is therefore that the higher-order frequency 
derivatives we observed are manifestations of a glitch recovery. 
In the AXP IRXS J 1708-4009, for instance, negative second fre- 
quency derivatives in the -(0.01^1. 3) x 10~ ^° Hzs~^ range has 



been detected jus t after a glitches jPib. Kaspi & Gavriill BoO^ : 
llsrael et al]|2007b|) . We found no evidence for glitches in the time 
spanned by our observations, but this does not exclude the possi- 
bility that such events occurred before the first observation of the 
outburst. 

Another possibility is that the frequency second derivative we 
measured (or at least a significant fraction of it) is linked to the 
magnetospheric activity of SGR 1833-0832. In the magnetar sce- 
nario the spin derivative is in fact expected to increase while the 
magnetospheric twist is growing, i.e. for instance in periods pre- 
ceding large outbursts, and to decrease in the aftermath (see e.g. 
Thompson, Lyutikov & Kulk arni 20 02; Bel oboro dov 2009). On the 
other hand, as pointed out bv lBeloborodovl ( f2009l) , a negative i> fol- 
lowing a period of bursting activity (as that detected here) can be 
still accounted for if a magnetospheric twist is suddenly implanted 
but its strength is moderate (twist angle less than ~1 rad). Then the 
twist may still grow for a while in spite of the luminosity released 
by dissipation monotonically decreases. Only after the twist angle 
has reached its maximum value, it will start to decay, together with 
the torque. 

Nevertheless, the timing properties of the source appear con- 
sistent with those typically observed in the AXP/SGR class. The 
spectral characteristics of SGR 1833-0832 are instead somewhat 
less usual. The X-ray spectrum is well described by a single black- 
body with temperature kT ~ 1.2 keV, definitely higher than 
what generally observed both in transient and 'persistent' mag- 
netars (kT < 0.7 keV, with the ex ception of SGR 0418-1-5729; 
lEsposito et al.l2010j ; lR"ea et alj2010h . We cannot however exclude 
that the high temperature of SGR 1833-0832 partially results from 
a bias in the spectral fits that, given the paucity of counts below 
~2-3 keV (owing to the large absorption), acts in the direction of 
increasing the temperature, in order to account for a second spectral 
component which cannot be properly modelled (see Section l372] ). 

The large increase in flux, a factor >20 above the quiescent 
level, is indicative of a rather powerful event, not dissimilar from 
those observed in other transient sources. This, however, seems to 
be in contrast with the low level of activity s een in SGR 1833-08 32, 
from which only few bursts were detected jGogtis et aLlboiOj) . A 
scenario in which we are presently witnessing the later stages of an 
outburst which was caught during its decay rather than at its onset, 
could account for the low activity but seems difficult to reconcile 
with the current high value of the temperature (if the thermal com- 
ponent in the spectrum of SGR 1833-0832 is indeed at kT ~ 1.2 
keV). A further possibility might be that SGR 1833-0832 experi- 
enced a (crustal) heating episode which, however, was not accom- 
panied by (or did not trigger) a twisting of the external field, as 
possibly suggested by the non-detection of a hard X-ray tail. In this 



iHobbs. Lvne & Krameill2010l) . Moreover, considering the correlation be- 
tween timing noise and spin-down r ate, one would expect for SG R 1833- 
0832 Ag to be approximately -0.6 jHobbs. Lvne & KrameJbOlOh . 
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case, however, an alternative explanation, not related to the twist 
evolution, for the presence of a negative v must be sought (perhaps 
a long-term postglitch recovery, see above). 

Immediately after the BAT trigger that led to the discov- 
ery of SGR 1833-0832, its observed flux was of ~3.8 x 10"^^ 
erg cm"'^ s~^, and it decreased by ~15% during the following 5 
months (which is not unusual for magnetars; e.g. iRea & Espositol 
l201lh . to X 10"^^ ergcm~^s"\ Although present data 

do not allow us to discriminate among different decay patterns, 
the flux evolution can be satisfactorily described by a broken 
power law, similarly to the case of SGR 0418-1-5729. The flux 
of SGR 1833-0832 remained fairly constant for about 20 days 
(ai ~ 0), when the decay index changed to Q2 ~ —0.5. Interest- 
ingly, alsoJhedecayofSGR 04 18-1-5729 exhibited a break at ~20 
days ( lEsDositoetal]|2()I0ah. As in the case of SGR 04 18-1-5729 
(see the discussion in lEsposito et al. I l2010al) . the presence of a 
break, together with the rather flat decay indices, might be diffi- 
cult to r econcile with the predictions of the de ep crustal heating 
scenario l lLvubarskv. Eichler & Thompsonl2002h . F oc with 
n ~ 2-3. This model pr ovides a satisfactory descr iption of the flux 
decay in SGR 1627^1 ( iKouveUotou et alj|2003 . although subse- 
quent analyses found no eviden ce for the plateau aroun d ~ 400- 
800 d predict ed by the mode l; iMereghetti et all liooel) and AXP 
IE 2259-1-586 fchu et al.ll2008l) . At the same time, however, it can- 
not explain the variety of decay patter ns which emerg ed from re- 
cent observations of SGRs/AXPs (e.g. lRea et ai]|2009l) . This may 
reflect the intrinsic limitations of the model (which relies an a sim- 
plified, one-dimensional treatment), or point towards the existence 
of different heating mechanisms which are at work during the out- 
bursts of magnetars, like in the case of SGR 1833-0832. 

Prompted by previous detections of a transient (pulsed) ra- 
dio emission following X-ray transient activity in other magne- 
tars jCamilo et al]l2006l |2007| ; iBurgav et al.ll2009h . we observed 
SGR 1 833-0832 with the ATCA and the 64-m Parkes radio tele- 
scope l lBurgav et al]|20IOl : Iwieringa et ai]|2010l) . No evidence for 
radio emission was found, down to flux densities of 0.9 mjy and 
0.09 mJy for the continuum and pulsed emissions, respectively. 
Similar upper limits on the pulsed radio emission, 0.1 mJy at 1.38 
GHz and 0.2 mJy at 2.28 GHz for a 10% duty cycle, were obtained 
with the Westerbork Radio Synthesis Telesco pe in the same day s 
(but not simultaneously with our observations : lG6giis et al .l2010ah . 
Despite this rather intensive coverage, the negative results of radio 
searches performed so far on SGR 1833-0832 cannot anyway be 
taken as conclusive because of the rapid variability of the pulsed 
flux shown by the known radio magnetars (see e.g. iBur gay et al.l 
l2009l) . Moreover, for a distance d = 10 kpc these limits translates 
into a pseudo-luminosity L — Sd? ^ 10 mJy kpc^, which is sig- 
nificantly smaller than the 1 .4 GHz luminosity of the other known 
radio magnetars at their peak (~ 100-400 mJy kpc^ : ICamilo et all 
[2006. 2 007. ; .Levin et al..20ld) but still much larger than the lumi- 
nosities of some known ordinary radio pulsarsl^H 



6.2 PSR B1830-08 and HESS J1834-08 

The discovery of the X-ray PWN generated by PSR B 1830-08, 
suggestive of a pulsar bow-shock and a possible diffuse emis- 
sion trailing the pulsar proper motion, provides new elements to 



See the online version of the Australia Telesco p e Nat ional 
Facility (ATNF) pulsar catalogue ^Manchester et al.l l2005h at 
|http://www.atnf.csiro.au/research/puIsar/psrcat/] 



the identification of HESS J1834-087. Accounting for a system- 
atic 20% uncertainty, the distance of P SR B 1830-08 is c ompatible 
with the one of W41, 4.0 ± 0.2 kpc frian et al.|[2007h . which is 
likely associated to HESS J 1834-087 given the low probability of a 
chance superposition. The pulsar spin-down age is also consistent 
with the dynamical age of remnant, roughly 1 00 kyr dTian et al] 
l2007h . As pointed out bv lMukheriee et al.N2009h , the pulsar moves 
fast enough to have reached its current position starting only 40 
kyr from the geometric centre of the remnant (see also Fig. |6] left 
panel). 

However, the Gamma-ray to X-ray flux ratio, 
F(E > 1 xeV)/-F'(2-io koV) ~ 50 for the compact emission 
and ~10 for the total X-ray emission, is lower than the one 
expecte d from the inverse correlation between this quantity 
and E jMattana et alj l200Si) , which should be around 600 for 
PSRB1830-08. If HESS J 1834-08 is the relic PWN generated by 
PSR B 1830-08, it is not very efficient in radiating in TeV gamma 
rays. Otherwise, it may be unrelated to PSR B 1830-08. The non- 
detection of a more extended X-ray emission trailing the pulsar 
does not support the association with HESS J 1834-087 either. 
Higher resolution and longer X-ray observations are needed to 
establish the nature of the halo and the detailed morphology of the 
compact emission, as well as singling out a possible contribution 
from the pulsar magnetosphere to the compact emission. 
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